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Abstract 

The method of exhaustion is generahzed to a simple formula that can 
be used to integrate functions under very general conditions, provided that 
the integral exists. Both a geometric proof (following the usual procedure 
for the method of exhaustion) and an independent algebraic proof are 
provided. Applications provided as examples include use of the formula to 
generate new series for common functions as well as computing the group 
velocity distribution resulting from waves diffracted from an aperture. 



1 Introduction 

The method of exhaustion was used by the ancient Greeks to integrate a limited 
number of functions before the development of calculus [1]. It will be shown here 
that this ancient method can be generalized to a simple formula to integrate 
functions under very general conditions, assuming that the integral exists over 
the interval in question. The formula is 



f.b oo 2^ — 1 

f(x)dx = {b-a)Y,Yl (-l)'"+'2-"/(a + ™(6 - a)/T) (1) 

^ n—l m—1 



The proof is simple and geometric in nature, directly following the procedure 
originally used for the method of exhaustion. An independent algebraic proof 
is also provided for the special case for which a = 0. 

2 Geometric Proof 

The geometric proof is a generalization of the successive approximation method 
historically used to obtain the exact area under the simple curves previously 



integrated by this approach. Consider an arbitrary function f(x) that is piece- 
wise continuous on [a,b]. The procedure to find the area under f(x) on [a,b] 
involves successive approximations with triangles that first intersect f(x) at the 
midpoint, x = a + (b - a)/2, and then at the new midpoints, x = a + (b - a)/4, 
X = a + 3(b - a)/4, and so on. The first such approximation is shown in figure 
1. Its area is 



Ai = ^{b-a)f{a+{b-a)/2). (2) 

The second such approximation adds two new triangles (figure 2), each shar- 
ing two vertices with the original triangle, and each bisecting one of the remain- 
ing arcs of f(x) at its third vertex. The two triangles have the area 

A^ = l{b-a) {/(a + (6 - a)/4) - /(a + (6 - a)/2) + /(a + 3(6 - a)/4)} . (3) 



A third such approximation (figure 3) adds the area 

^3 = l{b-a){f{a+{b-a)/8)-fia+{b-a)/A) (4) 

+/(a + 3(6 - a)/8) - f{a +{b- a)/2) + f{a + 5(6 - a)/8) 
-/(a + 3(6 - a)/4) + /(a + 7(6 - a)/8)}. 

Each new iteration n leads to 2"-l new triangles adding area A„ to further 
refine the area estimate. This procedure is continued indefinitely to "exhaust" 
the remaining area, leading to the formula 

„5 oo oo 2^ — 1 

/ f{x)dx = Y,An = {b-a)Y,Yl (-l)'"+'2-"/(a + m(6 - a)/2"). (5) 

•'^ n=l n=l m = l 

The special case for (5) when a = is as follows: 

,.b oo 2"-l 

/ /(x)dx = 6^ ^(-ir+i2-"/(m6/2"). (6) 

•"^ Ti=l m=l 

The method of exhaustion will converge to the value of the integral at least 
as fast as a geometric series, because when each new triangle is small enough so 
that the local curvature between intersection points on f(x) is slowly varying, 
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it will exhaust more than 1/2 of the remaining area left from the previous 
approximation. 

There are an infinite number of possible variations on this procedure. For 

example, each new iteration may involve a different number of triangles to ex- 
haust the remaining area, or each new triangle could intersect the function at a 
location other than the midpoint. 

3 Algebraic Proof 

An independent proof for (6) can be obtained starting with the identity [2] 
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This expression can be rewritten as follows; 
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or 
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Since 




(10) 



it follows that 




(11) 



Defining 




(12) 



3 



and noting that 



fc=o 



n=l m=l fc=0 

X cos{7nha/2'^) 

oo 2^ — 1 oo 
n— 1 m— 1 fe=0 

X sin(mfoa/2"), (13) 



it follows that 



Setting 



leads to 



oc 2"-l 

bJ2Y. (-l)'"+^2-" [G(m6/2") + G{-mb/2"')] 

n=l m— 1 

X cos(m6a/2") 

oo 2"-l 

~^^E E (-l)"+^2-" [G(m6/2") - G(-m6/2")] 

n— 1 m— 1 

X sin(m&a/2"). (14) 



G(a;) = H{x)e'' 



(15) 



/5 cxD 2" — 1 
H{x)dx = & 51 E (-l)™+^2-" [iJ (m6/2") + i?(-m6/2")] . (16) 
-6 „_i ^ 1 
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Note that H(x) can be in the form of a Fourier series which sums to the 
following: 



H{x) 



0: -b <x <0 
f{x);0<x<b /' 



(17) 



so that 



oo 2"-l 



/ /(x)rfa; = 6^ ^(-l)'"+i2-"/(m6/2"). (18) 

•^0 n=l m=l 



4 Applications 



The most immediate applications are series expressions for common functions 
that might otherwise be very difficult to derive. Some examples are presented 
below: 



sinx = a; ^ ^ (_l)™+i2-» cos(mx/2") ; (19) 

n=l m=l 



oo 2^-1 

cosx 

n=l m=l 



= 1 - a; ^ ^ (-l)"+i2-" sin(ma;/2") ; (20) 



/' '^dx = ± sin(m6a/2") ; (21) 



n=l m=l 



oo 2'^-! 



= 1 + X ^ ^ ^-l)^+-^2-''e'^'^/'^^ ; (22) 



n=l m=l 



oc 2"-l 

i(m6)^/4'^ . 



/ e-"^'dx = X! (-l)™+i2-"e-"( 

"'O n=l m=l 



oo 2"-l 



(23) 



Ina; = V V r 7V^; (24 

n— 1 m— 1 ^ ' 

(0 < X < oo) 



p\ = [ln(l/x)]^dx = ^ ^(-l)"+i2-"[ln(2"/m)]''. (25) 

•^0 n=l m=l 

(0 < p< oo) 

In addition to the above expressions, this method can lead to new insights 
into certain physical problems. One example involves the diffraction of waves 
though a two-dimensional aperture in an infinite screen. This problem has been 
solved exactly in integral form [3-5] so that the field at every point depends on 
the Fourier transform F{'kx,ky) of the aperture at z=0: 
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(26) 



Note that (26) satisfies the Helmholtz equation, 



V^(?i + k^(t> = (27) 

everywhere (k = Wo/c) as well as the boundary conditions on the screen. 

Under suitable conditions, the extension of (5) to convergent improper inte- 
grals can be made as follows, i.e., 



/ /(:c)dx = 6^ ^ ^(-l)'»+i2-"/(p6 + m6/2"). (28) 

''^ n=l m=l p=0 

Equation (28) is a result of breaking the integral into a series of definite 
integrals on [pb,(p-|-l)b]. The result is valid when each integral 



r(p+i)6 

exists and when 



/■(P+1)6 

/ f{x)dx 

Jvb 



[ f{x)dx (29) 
Jo 

tends to a finite limit L as B— !■ oo. 

With regard to (26), however, it can be noted that propagating waves will 
only occur when the transverse wavenumber (k^+k^)^/^ is lower than the cutoff 
wavenumber, so in the far field the integral can be evaluated to ±k with good 
accuracy as follows (assuming that F(k3;,kj^) is an even function with respect to 
both kx and k^): 



„,2 oo 2"-l oo 2"-! 

cPix, y, z, t) ~ — e-'""* E E E E (-l)"+«2-"-PF (mfc/2", qk/2P) 

n— 1 m— 1 p— 1 q—1 

X cos(m:rfc/2") cos(g2/A:/2f)e'"'=Vi-™V4"-9V4'' (30) 
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It can be clearly seen from (30) that the resulting field is due to a summation 
of an infinite number of plane waves, each propagating at a different phase 
velocity (and hence a different group velocity) based on its value of k^; and 
k,^. Thus, the process of diffraction leads to a continuous distribution of group 
velocities, having an amplitude distribution governed by the Fourier transform 
of the aperture. This insight becomes particularly clear upon application of (5) 
to expand (26). Such a result has an effect on Dopplcr shifts both in acoustics 
and in electromagnetic wave propagation, effects that are not otherwise apparent 
without the use of the generalized method of exhaustion. 
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